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Abstract 

We present an extended discussion of the previously noted possibility to extract 
\Vub\ from an analysis of the hadronic recoil mass spectrum in i? —»• Xu£u decays. 
Invariant mass spectra containing perturbative as well as nonperturbative corrections 
are given; their shape is manifestly sensitive to the three basic quantities pk 
and as, whereas the total integrated rate is much less so. Only a small fraction of 
b ^ u transitions generates a recoil mass Mx of at least M^. Moreover we hnd 
that the fraction of events with Mx < 1.5 GeV (to reject leakage from b ^ c due 
to measurement errors) exhibits fairly little dependence on mb and \Vub\ can 
then be extracted in a largely model-insensitive way. This conclusion is based on the 
applicability of the OPE to actual semileptonic B decays. A direct cross-check of 
this assumption and a determination of the required basic parameters of the heavy 
quark expansion will be possible in the future with more experimental data. 


Permanent address 




1 Introduction 


The KM parameters are fundamental quantities in the Standard Model. They have 
to be determined as reliably as possible for two reasons, one of a theoretical and 
one of a more phenomenological nature: (1) It is hoped that a future more complete 
theory will enable us to calculate these parameters. (2) CP asymmetries observable 
in strange and beauty decays are predicted in terms of these parameters. 

Numerical precision in the extraction of these quantities is highly desirable. The¬ 
oretical schemes involving quite different dynamical scenarios at very high energies 
tend to yield KM parameters that do not differ very much when probed at or below 
the electroweak scale. Also, since some of the predictions for CP asymmetries in 
B decays can be made with high parametric accuracy, one wants to translate this 
achievement into high numerical precision. Such considerations suggest a benchmark 
of better than 10% in accuracy for \Vch\ and \Vub\ as a desirable goal. 

With the emergence and increasing sophistication of heavy quark expansions (for 
a detailed evaluation, see |^) one has been able to translate ever more precise mea¬ 
surements into extractions of \Vcb\ with a theoretical uncertainty of about 5% and 
improving - something that could not have been expected a few years ago. 

The situation is much less satisfactory for \Vuh\- We know certainly that |14b| % 0 
holds since (i) the decays B —>■ p^v have been identihed and (ii) B —>■ X £u 

has been observed with lepton energies that are accessible only if X does not 
contain a charm hadron: Ei > (M^ — M^)/2Mb = 2.31 GeV. However to translate 
these hndings into reliable numbers concerning |Wb| is a much more difficult task 
theoretically: On the one hand the exclusive decays B —> n£i>, p£v depend on bound- 
state effects in an essential way (heavy quark symmetry can be relied upon here to a 
considerably lesser degree than in B ^ and although quark models have been 

employed, there is no reliable way for gauging their theoretical uncertainties). On the 
other hand in analyzing the endpoint spectrum for charged leptons in the inclusive 
decays B ^ X £v one encounters different sorts of systematic problems. Only a fairly 
small fraction of the charmless semileptonic decays B Xu£n, namely around 10% 
or so, produce a charged lepton with an energy beyond that possible for B Xc£n; 
in addition, the b ^ c rate is so much bigger than that for b ^ u that leakage from it 
due to measurement errors becomes a serious background problem; furthermore the 
endpoint region is particularly sensitive to nonperturbative dynamics. 

The importance of one such effect, namely the motion of the decaying heavy quark 
inside the hadron, was recognized a long time ago and the concept of ‘Fermi motion’ 
was introduced into phenomenological models, albeit in an ad-hoc fashion p, |^. It 
was later pointed out that Fermi motion emerges naturally in a dynamical treatment 
that is genuinely based on QCD implementing the heavy quark expansion through 
an operator product expansion (OPE) [Q. Conceptually it is similar to leading twist 
effects in deep inelastic scattering (DIS). Yet some subtle, though signihcant peculiar¬ 
ities arise making it somewhat different from the simple-minded treatment in quark 
models. In particular, it is not permissible to identify the basic expectation value 
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(to be defined below) with what is usually called the average Fermi momentum p^p, 
extract it from a fit to 6 —>■ c transitions in the most popular AC^M^ model, and then 
apply it at face value to 6 —>■ u decays. For p^p is in general different from even in 
the context of the AC^M^ model itself; furthermore QCD dictates using a somewhat 
different set of constraints than were historically used in the AC^M^ model [^. 

Nevertheless, we have learned more than just negative lessons. We know how 
to express the total width reliably in terms of \Vub\ and we have 

found descriptions of the Fermi motion that are - while not unique - at least fully 
consistent with everything we know about QCD. The theoretical tools involved have 
been described before p| . Here we want to concentrate on how they can be applied in 
a practical way. We will briefly discuss F(H —>■ before analyzing in detail how 

to extract \Vub\ from the hadronic recoil mass spectrum in semileptonic B decays. 


2 Total Semileptonic Widths 


The semileptonic widths F(H —>■ X iu) have been calculated through order l/mf 
0.1 i; the effect of the cubic terms is particularly small, and we neglect it in 
the rest of our discussion. The leading nonperturbative corrections are expressed 
through expectation values pQ and /i^ of the chromomagnetic and kinetic heavy 
quark operators, respectively 

l/a={B\b^aGb\B)/2MB, 1^ = , ( 1 ) 


with aG = where Gai 3 is the gluon field strength tensor; D denotes the covari¬ 

ant derivative. The measurement of Fsi(i?) has allowed the most reliable extraction 

of iWbl i: 
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The largest theoretical uncertainty resides in the value of pi as shown explicitly. 
For the total width depends mainly on the difference mb — me the error of which is 
controlled by pi] the remaining uncertainty in is stated in Eq.(g). The fourth 
term in the last bracket represents an estimate for the unknown effects of order l/rug 
(and higher) and of possible deviations from quark-hadron duality. The perturbative 
error reflects the uncertainty in the value of as and the weight of the higher order 
corrections. 

T{B ^ luXu) can be treated in complete analogy, and one finds Bi: 
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The dependence on /i^ has practically dropped out; the uncertainty introduced by 
rriiy is larger than in the b ^ c case, but still quite small. 

The real and highly nontrivial challenge is then of an experimental nature, namely 
how to measure r(i? —> The only feasible way would presumably be to find 

a kinematical discriminator between b ^ u and b ^ c transitions. This will be 
discussed next. Eq. (^) shows that the main uncertainty will be experimental rather 
than theoretical. 


3 The Hadronic Recoil Mass Spectrum 


As stated in the Introduction, phenomenological models include the motion of the 
heavy quark through a given ansatz. It was also recognized that a measurement 
of the hadronic recoil spectrum in semileptonic B decays 
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would offer intrinsic advantages for disentangling the b ^ u transitions relative to the 
conventional analysis of the charged lepton spectrum dT/dEf^{B —>■ Xlv) - if it could 
be done experimentally. In the lepton spectrum one can separate b ^ c and b ^ u 
kinematically only in the small slice 2.31 GeV < Ei < Mb/‘2, ~ 2.64 GeV, and not 
surprisingly ~ 90% of 6 —>■ m is hidden underneath the dominant b ^ c transition. In 
Mxi on the other hand, kinematical separation can be achieved in the much larger 
range < Mx < One might expect - and simple quark model computations 
like those of bear out the fact - that the bulk of the 6 —>■ m contribution lies below 
that for 6 —»• c when expressed in terms of Mx- Yet a more sophisticated analysis is 
required to see to which degree this holds true. The phenomenological descriptions 
involve some ad hoc assumptions of uncertain numerical reliability, and they suffer 
from obvious failures or at least limitations, for example in their description of the 
hadronic recoil spectrum for 6 —>■ c : the observed mass spectrum dominated by the 
two narrow D and D* peaks does not emerge from the quark model descriptions - 
instead these models yield very smooth functions for b —>■ u and b ^ c channels 
alike. Yet the QGD-based treatment yields a different expansion in the two cases 
rT| . The best tool for these studies is provided by the heavy quark expansion; first 
we will sketch this theoretical technology as it applies here and then analyze 6 —> m 
transitions in detail. 


3.1 The Methodology 

The theoretical treatment of inclusive semileptonic decays of B mesons carries a 
distinct similarity to deep inelastic lepton-nucleon scattering - this analogy can be 
pursued at great length. One defines a hadronic tensor as the B meson expectation 


3 











1 


(5) 


value for the transition operator 



= i / d^xe 4 ( 0 )} , 4 = qj^{l - 75)6 . ( 6 ) 
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The hadronic tensor can be decomposed into hve different Lorentz covariants 


h,,u{q) = -hg^^iqo, q^) + ^ 2 (^ 0 , q^)v^,v^ - ih{qo, q^)e^,yapVc,q 0 + 
+ 4 (^ 0 , q^)qiiqv + ^ 5 (^ 0 , q^){qi,Vu + quv^,) 


(7) 


{y is the 4-velocity of the decaying B), from which one obtains structure functions in 
the usual way: 


Wi{qo, q^) = 2Im4(go, q^) • 


( 8 ) 


All inclusive observables can be expressed in terms of these Wi. For Z = e or /i with 
m; ~ 0 only wi, W 2 and are actually relevant. 

Applying the OPE to the product of currents in Eq. one can express the 
structure functions through an inhnite series of expectation values of operators of 


higher and higher dimension. The coefficients become singular when one approaches 


free-quark kinematics. Therefore, even the limit rrib —>■ cxo does not allow one to 
evaluate the structure functions completely - it would require the resummation of 
inhnite series of equally important nonperturbative contributions. 

Nevertheless, with the large mass of the b quark one can pick up the leading 
operators for a given type of singularity, a procedure similar to resummation of the 
leading-twist contribution in DIS. Their effect is combined into the heavy quark 
distribution function which replaces the Fermi motion wavefunction of the quark 
models. 

For 6 —>■ c transitions (i.e. for g = c in Eq. (^) we have one more tool at our 
disposal: heavy quark symmetry imposes certain constraints on the properties of the 
structure functions. It can be revealed directly in the QCD l/mq expansion through 
studying the small velocity limit to derive SV sum rules without an a priori 
appeal to the underlying picture of strong interactions. Thus the tools are prepared 
to calculate (among other things) the hadronic recoil mass spectrum for b —>■ c in & 
way that is fully consistent with QCD and in agreement with the data - in contrast 
to what phenomenological models yield. One should keep in mind, though, that even 
the full power of the heavy quark symmetry still leaves large room for variations when 
the effects of higher order in l/rric and/or the recoil velocity, are addressed. 

With respect to 6 —>■ u the situation is a priori less favorable since heavy quark 
symmetry cannot be applied to the hnal state in 6 —>■ u and a small velocity treatment 
would make no sense! Nonetheless, in terms of the “global” characteristics relevant 
for inclusive decays, the heavy quark expansion yields many constraints. 

Similar to DIS, the moments of the structure functions, i.e. their integrals with 
integer powers of go at hxed g^, are given in terms of the expectation values of local 
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heavy quark operators of increasing dimension. De-convoluting these moments one 
can determine a heavy quark distribution function F{x) which encodes the effects of 
the Fermi motion. The dimensionless parameter x plays the role of the primordial 
momentum of the b quark normalized to A = Mb — rrib] x varies between —cx) and 1. 
All structure functions in the leading approximation can then be expressed in terms 
of F{x): 
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Here is a parton structure function including perturbative corrections. The 

exact form of this relation is not unique as long as only the leading-twist effects are 
summed up. Nevertheless, this particular form has some advantages, one of them 
being its transparent physical meaning (although it might not be quite obvious at 
hrst glance). The limitations - both of theoretical and practical nature - due to 
discarding the subleading-twist contributions and neglecting the actual dependence 
of F{x) on the velocity of the hnal-state hadronic system, are discussed in [^, . 

As stated above, the moments Oj of F{x) are given by the expectation values 
of local heavy quark operators. In practice we know only the size of the first few 
moments; in the adopted normalization one hnds 


Oq — 1 , Ui — 0 , 02 — 


/< 

3A2 


( 10 ) 


One then chooses a specihc functional form for F{x) and adjust its parameters to 
reproduce the phenomenologically deduced moments. This procedure was performed 
in Refs. |]1^ [T^ where the following ansatz was employed: 

F{x) = 6{1 — — x)°'[a + h{l — xY] . (11) 


With abundant experimental information one can in principle measure F{x), for 
example, in the inclusive decays 6 —>• s -|- 7 . 

Knowledge of the structure functions allows one to calculate all inclusive distri¬ 
butions. For example, the invariant mass squared of the hadronic hnal state has the 
following kinematic form: 

= Ml + q^ - 2MBqo = [ml + q^ - 2mbqo\ + 2{mb - qo)A A^ (12) 


Note the explicit term proportional to A; it manifests the l/rub nonperturbative effect 
in the average value of Mx 0- 

has a very intuitive quark model interpretation | 0 . 


Eq. 


Yet there are some 


relevant subtleties that have to be kept in mind for proper understanding as explained 
in detail in 


13, 14 
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• The nature of F{x) changes even qualitatively when going from heavy quark 

masses, > Arrib - to light ones - ^ A^. One of the advantages of the 

representation is that the relations stated in Eq. (pT]) for the first three 
moments hold for an arbitrary quark mass in the final state; the high moments 
essentially depend on it (as well as on through the quark velocity). 

• In the presence of gluon bremsstrahlung the hadronic recoil mass Mx for b ^ u 
can become much larger than \jKmb which would hold in the simplest quark 
picture. 

• We have adopted Wilson’s prescription for the OPE where an energy scale fi 
is introduced to separate long- and short-distance dynamics; the former are 
lumped into the matrix elements of the local operators while the latter are 
incorporated into their coefficients. No observable can depend on /r; yet as a 
practical matter it has to be chosen such that perturbative as well as nonpertur- 
bative corrections can be brought under theoretical control. More specifically, 
the moments Oj, the dimensionful parameter A as well as the distribution func¬ 
tion F{x) itself depend on /i. Some care has to be applied in keeping track of 
the /i dependence. 


While the second point is completely obvious, the first one is not and is actually quite 
mysterious from the perspective of a quark model description. 

Putting everything together we arrive at 
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A few more notes are in order on how the perturbative corrections are treated. We 
follow here Ref. [jl^ and include the exponentiated one-loop corrections at a given g^, 
evaluated with a fixed rather than a running ag, without an explicit infrared cutoff 
0 . The reasons behind this approach are given in |0. To be consistent, we then 


employ the b quark pole mass computed to one-loop accuracy with the same (frozen) 
coupling: 

, . 16 

mb ~ mR/i) -F --/i (14) 


9 TT 


and likewise for the kinetic term 


• (15) 

O 71 

In principle, the complete perturbative coefficient functions and the distribution func¬ 
tion separately do not allow considering the limit qi —>■ 0 whereas only the convolution 
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(p!3|) is /X-independent. Within the accuracy of our calculations, however, it does not 
pose problems. In particular, discarding the running of allows one to perform 
integration over gluon momenta down to zero. In respect to the moments of F{x) we 
considered above - as pointed out in Ref. - this corresponds to subtracting the 


would-be perturbative contribution to the expectation values of local operators from 
the domain below /x. This contribution must be evaluated in the theory with frozen 
coupling. This subtraction is performed in the equations above. 

The approximations we made in the treatment of the perturbative corrections are 
not expected to produce signihcant distortion. Even in 6 —s -f 7 decays for the 
actual b quark mass the major impact was found to be due to the soft primordial 
distribution [jT^. Moreover, even including the effects of running the resulting pho¬ 
ton spectrum (which measures the recoil mass spectrum) changed only a little upon 
variation of /x and upon incorporating the running of ag - while short-distance and 
long-distance parts separately changed radically. The stability in the case of semilep- 
tonic decays must be even better, since due to the sizable average value of q‘^, the 
effective energy release is much smaller and logarithms of the ratios of the momentum 
scales encountered in the perturbative calculations become rather insignihcant. 

To summarize our brief review of the relevant methodology: 


• The observable distributions are described in terms of structure functions which 
in turn are expressed - to leading approximation - through the convolution of 
a short distance rate and the heavy quark distribution function F{x). 

• F{x) encodes the main impact of nonperturbative dynamics conventionally re¬ 
ferred to as Fermi motion. In principle, it can be reconstructed from its mo¬ 
ments. In practice, only the hrst few moments are known. One then chooses 
some reasonable ansatz for F{x)-, its parameters are adjusted so as to repro¬ 
duce its known moments. While one cannot claim that this specific form for 
F{x) is derived from QCD in a unique manner, it represents a dynamical and 
self-consistent realization of QCD and its known constraints. Once more con¬ 
straints become available, they can in turn be implemented; i.e., the ansatz can 
be rehned step by step. 

• While many results or expectations previously inferred from phenomenological 
models re-emerge, it would be quite inappropriate to say they were ‘reproduced’. 
Deriving them from QCD proper represents significant conceptual progress; it 
has also warned us against various potential pitfalls and sharpened our vision 
regarding various subtleties that are quite signihcant even quantitatively, yet 
had escaped notice before. 

• The shape of the b ^ u distributions in practical implementation of this pro¬ 
gram is controlled by three quantities, namely 

— mb, the 6 quark mass ; 

— /x^, the kinetic operator for the b quark; 
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the strong coupling, Og- 


3.2 Analysis of 6 ^ w Recoil Mass Spectrum 

The tools appearing above are now applied. We use as central values for the basic 
parameters 

m6 = 4.82GeV, /f^ = 0.4GeV% a, = 0.3 (16) 

corresponding to the parameters of F{x) 


[a, b/a, c] = [0.5, 0.39, 1.75] 


(17) 


in Eq. (pdj) ; k = 1 has been preset. The resulting shape of the hadronic recoil 
spectrum is shown in the solid curve in Fig. 1. Adopting = 0.1 and = 0.5 instead 
produces the dashed and dotted curves which do not differ radically. The fact that 
dr /dMx shows such low sensitivity to the value of provides us with an a posteriori 
justihcation for our usage of a hxed coupling. This fact was expected. The main effect 
of the running of can come from the low-scale part of the bremsstrahlung. The 
exponentiation of the soft corrections, however, suppresses the domain where they 
are most pronounced. 

With hxed ag the effect of the exponentiation of the soft contributions is not very 
strong. For example, in the purely perturbative calculation of the hadronic mass 
distributions it shifts the value Mx, below which 60% of decay events are expected, 
from 0.55 GeV down to 0.4 GeV. This effect of softening the M^-distribution must 
be more prominent if one uses larger or running ag. 

As anticipated, the mass spectrum is very broad and extends even beyond Mo 
- yet only a small fraction does so, namely ~ 10%. Due to measurement errors 
there will be a tail from 6 —> c transitions below Mj^. To avoid this leakage one can 
concentrate on recoil masses below a certain value Mmax < Mu. The actual choice 
of Mmax is driven by competing considerations: the lower Mma x, the less leakage 
from 6 —> c will occur - yet the smaller the relevant statistics, as expressed through 
<l)(Minax), the fraction oib ^ u events with Mx below 
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There is a third consideration to be elaborated now. Since the recoil spectrum 
is shaped by nonperturbative dynamics, it is sensitive to the values of mb and /i^. 
Varying fil between 0.2 and 0.6 GeV^ while keeping ag = 0.3 and fhb = 4.82 GeV 
hxed one obtains the three curves of Figs. 3 and 4. Ghanging fhb between 4.77 
and 4.87 GeV, one obtains a somewhat larger difference in the spectrum - as shown 
in Figs. 5 and 6. With abundant statistics one can, in principle, distinguish the 
different curves and try to extract mb and - in addition to Yet it will take 

some time before such statistics become available. Realistically, we expect here only 




a cross-check of the consistency with other, dedicated experimental evalnations of 
these parameters. 

In the meantime one can qnite prohtably concentrate on analyzing $(Mmax); with 
the primary goal placed on determining iKifel- Since the theoretical nncertainty in 
the total width Ts\{b —>• u) or in the ratio rsi(6 —>■ u)/Ts\{b —> c) is small and depends 
on different type of effects, we dwell now on the nncertainty of predicting the rate 
in the kinematics where b ^ u decays can be experimentally disentangled from the 
KM-allowed transitions. In other words, we mnst nnderstand how accnrately one can 
calcnlate theoretically <h(Mmax) for realistic valnes of the cntoff mass Mma x. The hrst 
step is clearly to analyze its dependence on not yet precisely known parameters of 
the heavy qnark expansion. 

In Fig. 4 the predictions are plotted as a fnnction of M^ax for the three val¬ 
nes jll- = 0.2, 0.4, 0.6 GeV^ with fhb = 4.82 GeV hxed and in Fig. 6 for fhh = 
4.78 GeV, 4.82 GeV and 4.87 GeV with = 0.4 GeV^ hxed. We see immediately 
that the lower is chosen, the higher the sensitivity to long-distance dynamics 

becomes. This is not snrprising qnalitatively since for a lower Mm a x a smaller fraction 
of the overall rate is inclnded. 

The intervals of variation for rrib and we consider are qnite conservative. The 
actnal existing nncertainty in the latter seems to be 1.5-2 times smaller |^; it is hardly 
possible that the b qnark mass entering onr calcnlations can vary in a wider interval. 
Taking for orientation M^ax = 1-5 to 1.6 GeV we observe that the dependence of 
the fraction of the events $ (Mmax) on rrih is more prononnced. To illnstrate this 
dependence more clearly, we give the exploded version of Figs. 5 and 6 in Figs. 7 and 
8 where we vary mb by twice larger amonnt ±100 MeV. The variation in the fraction 
of the events remains reasonably small for Mm ax ^1.6 GeV. 

The stability deteriorates when one descends already a few hnndred MeV below. 
One not only loses statistics, the fraction $ becomes essentially dependent on the 
nonpertnrbative parameters. 

Thns we see that the goals of larger statistics and smaller sensitivity to nonper¬ 
tnrbative effects favor the selection of as high a valne for Mm ax as experimentally 
feasible, considering leakage from b ^ c. 

To be more specihc: taking the plots literally, abont 70% of all b ^ u events 
lie below Mm ax = 1.5 GeV with a spread of ±8%; for M^ax = 1-3 GeV and 1.0 GeV 
these nnmbers read ~ 60% and ~ 35%, respectively, with a spread of ±12% and 
±15%, respectively. These are nncertainties on top of those for the total B —> luXu 
width listed in Eq. (|^) which are enconntered in converting the semileptonic width 
into \Vub\- For M^ax ^ 1.5 GeV the nncertainties in the width pale in comparison to 
those dne to the shape of spectrnm. 

In reality, one mnst allow for additional nncertainties associated with other ap¬ 
proximations made in the analysis; most profonndly, the effect of the snbleading 
nonpertnrbative effects in the strnctnre fnnctions not captnred by naive Fermi mo¬ 
tion is qnite non-negligible for the actnal b qnark mass. Their theoretical discnssion 
and the estimate of related nncertainty can be fonnd in [IT^ M. As long as one does 
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not aim for an absolute precision in below ~ 10%, they are not expected 

to be important. Nonetheless, additional work toward better control of remaining 
uncertainties is clearly needed. Without a doubt, it will receive necessary theoret¬ 
ical attention as soon as the experimental feasibility of such measurements will be 
established and their ultimate precision is understood. 


4 Summary and Outlook 


We have shown here that if one succeeds in measuring the hadronic recoil mass spec¬ 
trum in semileptonic B decays, one can extract |Wb| in a theoretically clean and 
accurate way. In our theoretical analysis we split the problem into two parts - one 
is calculating the overall semileptonic 6 —> m width of B mesons, and the second is 
determining the fraction oi h ^ u events that can be discriminated kinematically 
against b ^ c processes. A good theoretical control over the hrst theoretical ingredi¬ 
ent allowed us to concentrate on the more subtle second aspect. 

Perturbative as well as nonperturbative QCD dynamics have a large impact on the 
mass spectrum - they broaden it considerably. Yet the numerical analysis suggests 
that they do not change very much one basic feature, namely that only a small tail 
distribution extends beyond Mo- Furthermore, while the shape of the spectrum is 
sensitive to long-distance dynamics, its integrated weight, say below 1.5 GeV is much 
less so. There is nothing magic about the value of 1.5 GeV. What is important 
is that one selects a cut-off that is high enough for most of the 6 —>■ m transitions 
to contribute below and low enough for almost all 6 —c leakage to occur above. 
Nevertheless, not only statistics, but also the issue of theoretical reliability calls for 
attempts to set as high a cut-off as possible. 

Thus the following scenario emerges: measuring the semileptonic decay rate below, 
say, Mmax = 1.5 GeV will enable us to extract |Wb| with a theoretical uncertainty 
that does not exceed 10 - 20 percent. An overall uncertainty of below 10% appears 
achievable or at least not impossible since with 10^ — 10® B decays and BR(i? —> 
IvXu) — 0.0015 one estimates that a statistical accuracy level ~ 0(1%) would be 
available, provided a large fraction of the decay events can be successfully utilized. 
As described above, a particular premium above and beyond statistics has to be 
placed on achieving a high value for Mmax, not below 1.5 GeV. 

Meanwhile, if the studies will demonstrate the feasibility of such experimental 
measurements, the theoretical description can be improved. The feasible method 
of obtaining accurate perturbative corrections to the decay structure functions was 
suggested in the so-called APS which combined advantages of exact one loop- 


expressions, exponentiation of soft infrared/collinear effects and incorporating of the 
running of the coupling. It intrinsically includes the normalization point and thus is 
free of problems associated with using the pole mass of b quark encountered in the 
“practical” OPE. Yet for the reasons elucidated above we do not anticipate qualitative 
changes in our predictions. In particular, the resummation of next-to-leading logs 
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does not seem very promising in the semileptonic decays due to a limited gap in the 
momentum scales. Without a cutoff, however, the resummation of leading logs seems 
technically necessary not to violate unitarity for soft emissions. It is important to 
check that the sensitivity to the exact treatment of the perturbative corrections is 
not high, since application of the perturbative expansion at a precise level for the 
soft gluons with ^ 1 GeV seems problematic. 

In the near future one will presumably determine rrih and /i^ from b ^ c transitions 
with even better accuracy than stated above. Lastly, with increasing statistics one 
can measure the b ^ u recoil mass spectrum with good accuracy, rather than its inte¬ 
gral as a function of M^ax- This in turn will allow us to independently place bounds 
on ruft and in a systematic way. Comparing these bounds with other determina¬ 
tions of these basic quantities will serve as an essential quality check on systematic 
uncertainties and the quantitative applicability of the heavy quark expansion. 

If accurate measurements of the recoil mass spectrum are feasible at all, it is 
advantageous to perform such studies separately for charged and neutral B. In the 
KM-suppressed semileptonic decays one can observe the effects of Weak Annihilation 
as the difference in the decay characteristics for charged and neutral B. Moreover, 
in the decays to massless leptons one is sensitive to the nonfactorizable parts of the 
expectation values of the four-quark operators responsible for the effects of Weak 
Annihilation |^. On the other hand, it is just the theoretical expectations for 


nonfactorizable parts which are a subject of some controversy in the literature at 
the moment. Gaining direct experimental information not only can put theoretical 
predictions of beauty lifetimes on more solid grounds, but also has an independent 
theoretical interest. Ultimately, Weak Annihilation can even introduce an element of 
theoretical uncertainty in the extraction of \ Vub\ in the discussed way unless its effect 
is controlled through such studies. 

At the same time, theoretically, the effects of Weak Annihilation are expected to 
come from small Mx (or even from small hadronic recoil energy E^) 0; singling 
out this domain can amplify the effect and make it detectable even if its overall scale 
is insignihcant. 
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Figure 2: The integrated fraction of the events 
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Figure 3: Dependence of dT/dM^ on /i^: long-dashed, solid, and short-dashed 
lines correspond to = 0.2 GeV^, 0.4 GeV^ and 0.6 GeV^. The 6-quark mass rhf, = 
4.82 GeV, as = 0.3. All distributions are normalized to the same total width. 



Figure 4: <F(MA for fll = 0.2 GeV^ 0.4 GeV^ and 0.6 GeVl 




Figure 5: Dependence of dT/dM^ on nib: long-dashed, solid, and short-dashed lines 
are fhb = 4.77 GeV, 4.82 GeV and 4.87 GeV. = 0.4 GeV^, = 0.3. 



Figure 6: <h(Mr) in the same setting as in Fig. 5. 




